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A macroscopic theory for the molecular or Casimir interaction of dielectric materials with ar- 
bitrarily shaped surfaces is developed. The interaction is generated by the quantum and thermal 
fluctuations of the electromagnetic field which depend on the dielectric function of the materials. 
Using a path integral approach for the electromagnetic gauge field, we derive an effective Gaussian 
action which can be used to compute the force between the objects. No assumptions about the 
independence of the shape and material dependent contributions to the interaction are made. In 
the limiting case of fiat surfaces our approach yields a simple and compact derivation of the Lifshitz 
theory for molecular forces For ideal metals with arbitrarily deformed surfaces the effective 
action can be calculated explicitly. For the general case of deformed dielectric materials the ap- 
plicability of perturbation theory and numerical techniques to the evaluation of the force from the 
effective action is discussed. 

PACS numbers: 12.20.-m, 33.90.-|-h, 42.50.Ct 

I. INTRODUCTION 

The last years have witnessed a resurgence of theoretical and experimental research on Casimir interactions between 
macroscopic and mesoscopic objects 0, 0, S IS IE 0- The most mentionable recent achievements include on the 
experimental side the high precision measurements of the attractive Casimir force between metallic surfaces |^ |^ 0, 
Ul, 12, 13 ■ ll^lTsllT^llTj and the simultaneous study of both the critical Casimir force due to (thermal) order parameter 
fluctuations and the electrodynamic Casimir force in superfluid Helium films near the critical point HaUJil^al- On the 
theoretical side there is an ongoing attempt to describe the Casimir interaction between non-ideal metals with finite 
conductivity, including simultaneously the effect of finite temperature, see [2l] | for a recent summary of the used and 
partly controversial approaches. Another direction of great interest is the study of the strong geometry dependence of 
the Casimir force which is inevitably linked to the non-additivity of fluctuation induced forces j22y2S, .24 .251] . Even a 
potential chan ge f rom attractive to repulsive forces due to either the material properties [2^ I2^l28l | or the geometry 
of the objects [29| has been discussed. As to the comparison between experiment and theory, Casimir's prediction 
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for the force Fq between ideal metallic and parallel flat surfaces of area A and distance H at zero temperature has been 
confirmed by the recent experiments within a few percent of accuracy. However, there is broad agreement that the 
experiments have shown that the inclusion of the material properties of the surfaces, of the roughness and geometry of 
the surfaces, and of finite temperature effects is indispensable. A solid theoretical account of these effects is necessary 
in view of the importance of the experimental results to the test of unified gauge theories of fundamental interactions 
[sif and the design of nanotechnological devices ,16, sj- 

Each of the modifications to the ideal prediction of Eq. |^ introduces at least one corresponding length scale in the 
interaction energy. If the interacting materials are not ideal metals, the electromagnetic fluctuations are not reflected 
perfectly at all wavelengths by the surfaces. Then one expects Eq. ^ to hold only for distances H much larger than 
the plasma wavelength Ap of the metal. In the recent experiments \p is of the order of O.l/xm. At smaller distances 
the force will be reduced to F = Tym-Fo by a material dependent factor ?7ni < 1 compared to the ideal force. Since 
the experiments are usually performed at room temperature, thermal fluctuations of the electromagnetic field tend 
to increase the force for separations H which are larger than the de Broglie wavelength of photons. At — fic/{kBT) 
(w 7/im at 300°K), leading to F = ?7t-Fo with the temperature dependent factor r/T > 1. Finally the presumably 
most difficult to calculate change of the force Fq comes from the surface geometry. The geometry modifications can 
be divided into two different types. The first is a in general unwanted stochastic surface roughness, the second an 
intentionally designed surface structure like a corrugation. The latter has been studied in a recent experiment to 
study the geometry dependence of the normal and lateral Casimir force Both types can be characterized by 

the deformation amplitude a and the roughness correlation or corrugation length A. Recent theoretical work for 
ideal metals at zero temperature has shown that the force is generally increased by deviations from the flat surface 
geometry, at least for uni-axial structures [2^ I2M |2^ . Therefore we can write F — rjgFo with a factor rjg accounting 
for geometry dependent changes of the force. 
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So far, we have discussed the different modifications of the ideal force Fq independently. One can expect that 
this is indeed justified if the characteristic length scales of the modifications are widely different from each other. 
However, this is by no means always the case in the recent experiments, especially at the point of closest approach of 
the interacting surfaces, and in nanotechnological devices. In view of this, it is of importance to consider correlations 
between the force modifications. That means that one has to assume a more general form of the force 

F^VmVT Ve (1 + Acorr)i^O (2) 

with a new term Acorr accounting for correlations between different modifications. The present state of the techniques 
available in literature, however, does not allow to determine Acorr in general. A commonly used but uncontrolled 
approximation is to set Acorr — 0. Most of the recent theoretical studies of correlation effects have been devoted to 
the simultaneous effect of finite conductivity of metallic surfaces and finite temperature, see for a recent overview 
and [s^l for a proposed experiment to measure these effects. Common to almost all of the existing theoretical work 
on the correlations between conductivity and temperature corrections is that it starts from the so-called Lifshitz 
theory for molecular interactions between macroscopic objects The Lifshitz theory provides a formula for the 
force in the rather general case of dielectric bodies at arbitrary temperature but with flat and parallel surfaces. The 
key difference in the approaches to correlation effects is the treatment of the zero Matsubara frequency term of the 
Lifshitz formula when it is applied to different models (Drude, ideal metal or free plasma model) for the dielectric 
function of the metals. In addition, there are a few approaches which are not based on the Lifshitz theory but employ 
a surface impedance boundary condition in order to account for the coupling between electromagnetic fluctuations 
and a metallic surface. To date, however, it appears that there is no broad agreement on a correct prescription for 
the evaluation of the Lifshitz formula in the case of non-ideal metallic surfaces at finite temperatures. 

The situation becomes even worse when one considers deviations from the parallel flat plate geometry and corre- 
lations between geometry, conductivity and temperature modifications of the Casimir force. To our knowledge there 
is no complete theory for correlations between geometry induced modifications of the force and the above discussed 
corrections available in literature. In fact, only recently an exact description for the sole geometry dependence of the 
force, making no additivity assumptions of two-body forces, has been given l23l l23|. So far, surface roughness 
in combination with finite conductivity has been studied only by neglecting correlations between both effects '35|. 
However, the characteristic length scales of surface roughness or designed surface corrugations in nanotechnological 
devices can be close to the relevant length scales of the material as, e.g., the plasma wavelength. Therefore it would 
be very useful to have an analog of the Lifshitz theory for more general geometries. Such a general theory should 
yield the Casimir interaction depending on the dielectric function of the material, temperature and a height profile 
describing the surface geometry as input parameters. The development of such a theory is the purpose of the present 
work. 

In this paper we will introduce a novel macroscopic approach to molecular forces between dielectric media. It will 
be based on a path integral technique for fluctuation induced forces which was previously developed for ideal metals 
at zero temperature 0, 12^ Is^ . Non-local boundary conditions for the electromagnetic gauge field are employed to 
treat the interaction between electromagnetic fluctuations and matter. The boundary conditions can be viewed as a 
reformulation of the so-called extinction theorem of classical electrodynamics '37*, '38', '39]. The important new property 
of our approach is that surface deformations can be included without any assumption about the correlations between 
contributions to the force from geometrical and material properties. We derive an effective Gaussian action which 
is a functional of the frequency dependent dielectric function of the material and the height profile of the surfaces. 
The effective action is a possible starting point for future detailed analyses of the effect of correlation on the force as 
described by the Acorr in Eq- (01 ■ We demonstrate the efficiency of our approach by looking at two particular limits 
of interest. First, we consider flat surfaces of dielectric media. In this case we obtain, as a byproduct of our theory, a 
compact and concise derivation of the Lifshitz formula for molecular forces Q in the language of quantum statistical 
mechanics. An even simpler derivation of this formula is found by a scalar field approach, cf. the Appendix, which 
should be compared to other derivations of the original result of Lifshitz. As second limiting case, we consider 

ideal metals with arbitrary deformations. Then the effective action assumes a simple form which can be determined 
explicitly. For general deformed dielectric media, the effective action can be used as a basis for perturbative [2^ 
or numerical |24| computations of the correlation term Acorr- 

The outline of the rest of the paper is as follows. In the following section we develop the path integral approach 
for deformed surfaces of materials which are characterized by a general frequency dependent dielectric function. We 
derive the non-local boundary conditions which describe the reflection and refraction properties of the interacting 
bodies. By integrating out the electromagnetic gauge field, an effective action for the interaction between the bodies 
is obtained. In Sec. IIIII we apply our theory to calculate the force between two flat surfaces of dielectric media. 
In this case, the known Lifshitz formula for molecular forces is found. The effective action for deformed surfaces of 
ideal metals is computed explicitly in section IIVI Section provides a discussion of the relevance of our results to 
perturbative and numerical analyses of the Casimir or molecular interaction between macroscopic objects. A rather 
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short and concise derivation of the standard Lifshitz theory in terms of a scalar field is left to the Appendix. 



We will develop a macroscopic theory which allows to calculate the interaction between materials of rather general 
shape. Instead of considering directly the field emitted by the fiuctuating dipoles in the material we view the interaction 
as occurring through the modifications of the quantum (and thermal) fluctuations of the electromagnetic field between 
the materials. No direct reference is made to the electromagnetic field fluctuations in the interior of the materials. 
The effect of the dipoles induced by the external fluctuating field will be described by material dependent boundary 
conditions which are defined on the surface of the material. Our method is based on a path integral quantization of 
the electromagnetic gauge field. This approach has full generality in the sense that it can be applied to any body, 
characterized by its dielectric function, with any surface profile, described by a height field, at any temperature. 

The common approaches for computing the force between materials is to first determine the solution of Maxwell's 
equations both inside and outside the materials, and than to evaluate the force either from the stress tensor or from 
the zero-point energy of the modes using the so-called argument theorem, see, e.g., Ref. The problem with these 
approaches is that they are not suited to treat non-flat surfaces since deformations in general lead to a complicated 
modification of the mode structure and make the solution to Maxwell's equations a hard task. In the following, we will 
formulate the interaction between deformed materials within the language of quantum statistical mechanics. Since 
this formulation makes no explicit use of the individual eigenfrequencies of the modes it is better targeted for the 
treatment of deformations. 

We consider the two interacting media as filling half spaces which are bounded by deformed surfaces S'a, a — 1,2. 
The deformations from a plan-parallel geometry of mean surface distance H are described by the height functions 
/iq(x||) with X|| the lateral surface coordinates, see Fig. ^ The media are characterized by their complex dielectric 
functions ea{uj), respectively. The gap between the media is assumed to be vacuum, e(uj) = 1. The free energy T of 
the photon gas in the gap between the two surfaces can be calculated from the imaginary time path integral for the 
electromagnetic gauge field A. In the absence of media, the vacuum partition function Zq is given by 



where we have introduced a complex valued gauge field which leads to a double counting of each degree of freedom. 
The reason for this will become clear below when we discuss the boundary conditions at the surfaces. The Euclidean 
action Se{A*A) is the imaginary time version of the action S{A*A) of the electromagnetic field in the Minkowski 



II. PATH-INTEGRAL FORMULATION OF MOLECULAR FORCES 




(3) 



A 




L - 



R 




£2(C0) 



H 







-L - 



FIG. 1: Two deformed surfaces Si and S2 of dielectric media with dielectric functions ei(aj) and e2(<i'), respectively, separated 
by a gap of mean size H along the a;3-direction. The meaning of the auxiliary surfaces Ri and R2 is explained in the text. 
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space time with coordinates X = {t,x) = {t,x\\,x^) , 

SiA:A) = -]^j^{F;^F^'^){X)--J^{d,A*^){d,A'^){X) (4) 

where the first term comes from the Lagrangian of the electromagnetic field F^^'^ = d^A'^ — d'^A^ and the second term 
results from the Fadeev-Popov gauge fixing procedure which assures that each physical field configuration is counted 
only once in the path integral over the gauge field. The parameter ^ allows to switch between different gauges; all 
gauge invariant quantities calculated from this action like, e.g., the Casimir force, do not depend on ^. In the following 
we will use the Feynman gauge corresponding to ^ = 1. The coefficients in the action of Eq. Q differ by a factor of 
1/2 from the conventional definition of the action for a real valued gauge field in order to obtain the correct photon 
propagator which in Feynman gauge reads G;.,^ = g^j^^/K^ with momentum K = (a',k) , if^ = K^K^ = — 
and Minkowski metric tensor g^^ = diag(l, — 1, — 1, — 1). The Euclidean action S-e{A*,A) is obtained from Eq. 
by applying a Wick rotation to imaginary time which amounts to the transformations t —iT, oj iC, and 
— > zA", A*^ — > iA*^ , while the remaining components A^^ remain unchanged [i^li^ . Since this transformation 
corresponds to the change 51^^ — > —5fj,i, for the metric tensor, the Euclidean action in momentum space becomes 

1 °° r 

Se{A*,A) - ^ E y^^*'^(Cn,k)^El^,(C„,k)A^(C„,k) (5) 



n— 



where we allowed for a finite temperature T by introducing bosonic Matsubara frequencies Cn = 27rn//? with /3 = l/T. 
The Euchdean Green function is given by QE.,^lu{C, k) = 5f^j,Q{C,, k) with Q{C,, k) = ((^ + k^)^^. 

In the presence of the two media the free energy of the photon gas in the gap between the media is obtained from 
a restricted partition function. The restrictions are due to the boundary conditions for the gauge field which are 
imposed by the dielectric properties of the media. It will be shown below that there are three boundary conditions on 
each surface Sa which we number by j = 1, 2, 3. Each of these conditions implies the vanishing of a non-local linear 
combination of derivatives of the components of the gauge field. Thus the mean distance i7-dependent restricted 
partition function Z{H) can be written as 



= 4 h{A*A) nn n 4/ q,{Cn;^.^)Ancn,^ 



(6) 



where we enforced the boundary conditions by inserting delta-functions for all positions x on (flat) auxiliary surfaces 
Ra which are placed at X3 = ±L with sufficiently large L so that the surfaces Sa are located between them, see Fig.^ 
The final result for the force between the media will be of course independent of L. The differential operators 
depend via both the dielectric function Cq, and the normal vector on the surface index a. Their actual form will 
be computed below. The interaction (Casimir) free energy per unit area of the two surfaces Sa is then given by 

HH)--^in[ZiH)Z^HH)], (7) 

where A is the surface area and (5 = 1/T the inverse temperature. Z^o is the asymptotic limit of Z for large H so that 
T is measured relative to two surfaces which are infinitely apart from each other. The force per unit area between 
the surfaces is then given by = —dT/dH. 



A. Boundary conditions 



In this section we will derive the boundary conditions at the surfaces of the dielectric media. The boundary 
conditions are based on the optical extinction theorem of Ewald 37] and Oseen [s^l, see also [s^- This theorem 
states that part of the electromagnetic field produced by the molecular dipoles inside a medium exactly cancels the 
incident field, while the remainder propagates according to Maxwell's equations in continuous media. Ewald and 
Oseen proved the theorem for crystalline media and amorphous, isotropic dielectrics, respectively, from a point of 
view of classical molecular optics. Later, Born and Wolf extended the theorem to more general classes of materials 'soj. 
A relationship between the extinction theorem and the Lifshitz theory of dispersion forces for continuous media has 
been pointed out by Milonni and Lerner . They use the fact that the extinction theorem permits a reduction of the 
multiple-scattering problem for the molecular dipoles to the solution of the wave equation for the gauge field A with 
appropriate boundary conditions. From this they conclude that the extinction theorem shows that the macroscopic 
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Lifshitz theory for continuous media correctly accounts for all multiple-scattering nonadditive contributions to the 
force between flat surfaces. We will demonstrate that these concepts are useful to describe the interaction of even 
deformed surfaces. 

In the following we use an (equivalent) reformulation of the extinction theorem as a non-local boundary condition 
which enforces the laws of reflection and refraction at the surfaces of the interacting media. Our derivation of the 
boundary conditions follows closely the approach outlined in |4It| | . Let us start with the Helmholtz wave equation for 
the magnetic fleld B inside a dielectric medium which occupies the volume V with surface S = dV, 



0. 



The propagation of the field inside the medium is described by the material Green function satisfying 

V^ + eiLu)LO^WiLu;^,^') = (5(3) (x-x'), 



(8) 



(9) 



where we allow for a general frequency dependent complex dielectric function e(w). Applying Green's theorem to the 
components of B and to Q"^, one easily obtains, using JSJ and 



x'es 



g^(a;;x,x')(n'- V')B(cc>,x') - B(cc>, x') (n' • V'^^(c^; x, x')) 



B(w,x) , xeV 
, x^V 



(10) 



with the normal unit vector fi of the surface pointing into the vacuum. We will make use of this result for the case 
where the position x is located on an auxiliary surface which is placed outside the medium so that the integral has 
to vanish. Using a succession of vector identities and the macroscopic Maxwell equations for continuous media, 
V-B = 0, VxB = —iuje{uj)'E, the condition that the integral vanishes can be transformed to 



x'es 



-iwe(w)(n' X E(w,x')) + (n' • B(w,x'))V' + (n' x B(tj,x')) x V' g''{uj;x,x') = 



(11) 



The usual continuity conditions for the electric and magnetic field at dielectric surfaces |46j without external surface 
charges and currents show that the three terms in Eq. have to be continuous across the surface. Therefore we 
can use the vanishing of the integral as a boundary condition for the electromagnetic field on the vacuum side of 
the surface. As a side remark we note that if we had started with the wave equation for the electric field instead of 
the magnetic field we had obtained a similar expression as that in Eq. However, due the discontinuity of the 

dielectric function across the surface the condition that the integral vanishes had not translated to the field on the 
vacuum side. Finally, in the case of an ideal metal, e(cj) — > oo, and the integral in Eq. Hill) is dominated by the first 
term. In this limit the integration can be carried out, leading to the well-known condition n x E = for ideal metals. 

The condition of Eq. l(TT)l can now be used to determine the differential operators (C, x, x') appearing in Eq. 0. 
We express the electric and magnetic field in terms of the gauge field. After a Wick rotation to imaginary time, the 
corresponding relations read Ej = —(A^ — idjA'^ and Bj = EjkidkA^ in Euclidean space. Multiplying Eq. (|11|) with 
(Ceq)"^ and decomposing /^"^(C; x, x') = n^(x')£*^^((^; x, x') with respect to the components of the normal vector 
by using the standard summation convention for fc, one gets after some algebra the explicit result 

/ 



£!"((; x,x') 



/:2"(C;x,x') = 



/:'"(C;x,x') 






--^[8382 - 8283] 


^[8183 + 8381] 


-^[8182 + 828,] \ 








C-^Jd282 - 8,8i] 


-182 




( + ^[8383-818,] 




-183 


-crAd3 


^8283 


C + ^[di8i- 8282]\ 





-^[8283 + 8382] 


-ck[^id3 -838i] 


^[9291-1-91 92] 


i8i 


C - Ck[^383 - 8282] 




ci^39i ) 


182 


crA92 


C - ^ [8i8i - 8383 


] -c7r^392 \ 


—i8i 


-( + ^[8282-8383 









^[8382 + 8283] 


-^[838,+8,83\ 


-^[828,-8,82]) 



g^i{C;x-x') (12) 



c;^°(C;x-x') (13) 



e^°(C;x-x'),(14) 



where fj^" is the Euclidean version of the Green function inside the medium which is defined by Eq. @ . The partial 
differential operators 8j are acting on the spatial argument of , whereas the "free" operators 8j are acting on 
the gauge field to which C"^ is applied. For non-deformed surfaces as described by the conventional Lifshitz theory, 
i.e., n = (0,0, ±1), only the last matrix is relevant. 
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B. General result for deformed surfaces 



Now we are in the position to calculate the partition function defined by Eq. © and by the operators in Eqs. I|12I14() . 
Similar to the approach of Refs. [3^l47| we introduce auxiliary fields in order to treat the delta-function constraints. 
However, here we will use complex valued auxiliary fields since the arguments of the delta-functions are complex 
in our problem. Moreover, the fields will be not defined on the original surfaces Sa itself but on the flat auxiliary 
surfaces i?^ since these are the regions on which the "external" positions x of the boundary conditions are located, 
cf. Eq. Introducing on each of the two surfaces Ra at X3 = La = (— with lateral coordinates xy the three 
fields ipajiCj X|| ), (j = 1, 2, 3), the delta-functions for fixed a and j can be written as 



nn4/ ^"m(c«;x.x')a'^(c,x') 



^ [^aj^aj] exp 



Xii Jx'£Sa 



{ raj (Cn , X|| )£'^^ (C„ ; (X|| , i„) , X') (C„ , X') + C.C. } 



(15) 



Inserting this representation in the partition function of Eq. ((^l, the complex gauge field can be integrated out, 
using the free action Se{A*,A) of Eq. jSJ. The partition function can then be expressed in terms of an effective 
quadratic action for the auxiliary fields, 



(16) 



with 



^eff[C,V'a,] =E / / C,(Cn,X||)M"'^'^'(C„,X||;C„',xI|)V;3i(Cn',xp, 

„ „/ Xii •/ Xm 



(17) 



where the usual summation convention applies to all indices. Since in Feynman gauge the propagator of A'^ is diagonal 
in fi, the resulting matrix kernel can be simply written as 



A^"'^'-'"'(C,X|i;C',x9 = 2^<5(C-C') 



5'^(C;(x|i,LJ,y) 4f (C';(xl,,i^),y') C?(C; 



y - y j, 



(18) 



where a summation over fi is implicit and ^(Ciy) is the free photon propagator with Fourier transform CJ(C,k) = 
(C^ -|-k^)~^. To simplify this result, and to prove the independence of the free energy on the choice of L, it is useful to 
rewrite the operators of Eqs. H12I14|I as C^°'{(^\ x, y) = ^^"'{Q G'e (C; x — y) in order to make their proportionality 
to the material Green function explicit. It is important to keep in mind that the differential operators ^'"^ act on 
the spatial arguments of Q'^^ as well as on those of the free propagator Q. Now the kernel in Eq. H18|) acquires the 
form 



Jy&SaJy'eSfi 

X (C; X - y)|,3.i„ (C; x' - y')|4=L, ^^(C; y - y'), 



(19) 



where ^'^^ acts on the primed coordinates and a summation over k and s is implicit. In momentum space, using 
t/|;(C,k) — {C'^e{iC) + k^)~^, the partially Fourier transformed material Green function can be written as 



e^°(C;k||,^) 



g-pa(C,k||)|z| 

2pa(C,kj|) 



(20) 



with PaiC, k|[) = y eQ(*C)C^ + k|. With this representation the kernel becomes 

g-p„(C,k||)|L„-i/3| g-P/3(C'.kJ|)|i/3-a3l 



M«^'^'(C,k|,;C',k') = 27r5(C-C') 



-ik,|y||-|-ik||y,| 



yes^Jyes, 2p„(C,k||) 2p4C',kj|) 

c n« hf [^^^ (C, k„ ) • ^'t«^ (C , kj, )] ^, c;(C; y - y') , 



(21) 
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where the differential operators ^'^"(C, ky) are obtained from the ^^"^ of Eqs. (|12I14(I by the replacements V|| = 

d-i) ik|| , ^3 (— 1)"Pq,. Thus the operators ^'^"(C, ky ) are acting via the remaining derivatives dj only on the 
spatial coordinates of the vacuum Green function Q(C^\y). At this stage it will become obvious that the free energy 
or force is independent of the positions 2:3 = ±i of the auxiliary surfaces i?Q,. Due to the construction of the surfaces 
i?Q,, we have |Lq, — yz\ = — 2/3)- A consequence of this is the important observation that the kernel can 

be factorized into 

A^"'''^"'(C,k||;C',kj|) = ,7„(C,k||)A^"'5.^'(C,k,|;C',kI|),7MC',kI|) (22) 
with the functions ria{C,i^\\) = cxp(— k|| )L)/2pq(C, k|| ) and the simplified L-independent kernel 

7W"''"''(C,k||;C',kj|) =27r^(C-C') / / e~*''ii''"+'''Ii^'i e"''"^^"^"^^'''"'^^^^"^^''*'''^^'''''!^^^] 

Jyes^ Jy'eSp (23) 

X nf [^'="(C, k|| ) . kj|)]^, g{(:- y - y'). 

From Eq. (|16|1 follows that the partition function Z{H) — det^^^^ M. with the determinant taken with respect to both 
the continuous (C, ky) and the discrete (a, j) arguments. Due to the structure of Eq. one has det oc det A^. 
Since the functions 77q,(C, ky) are independent of the mean surface distance H, the proportionality constant of the 
two determinants is independent of H, too. Therefore, this constant, and the L-dependence, will drop out of the free 
energy of Eq. |(7J) which can now be written as 

HH)^i^P^^'^'^^{mM-J), (24) 

where A^oo denotes ^A in the limit of asymptotically large H. The force per unit area between the two surfaces can 
be directly obtained from the kernel M. by 

F^-^TT(M-'dHM) (25) 

without the need to subtract the asymptotic expansion for large H. Here the trace has to be taken with respect to 
the Matsubara frequencies Cn, the lateral momenta ky, and the discrete indices a, j. Eqs. 124|) and (|25|l together 
with Eq. (|23|l represents the main result of our general approach. We will apply this formula below to specific 
model situations. Before proceeding along these lines, it might be interesting to discuss the above result. During 
the derivation of the result we worked within the Feynman gauge. This, however, poses no problem since the 
restricted partition function Z{H) can be considered as the expectation value of the boundary condition enforcing 
delta-functions with respect to the free action of the gauge field. The arguments of the delta-functions are composed 
of the electromagnetic field components, and are thus manifestly gauge invariant. Therefore, the kernel M must be 
gauge invariant. 

Let us first discuss the most simple situation where the kernel M is diagonal in momentum space so that the 
force can be calculated exactly. This will be the case when the geometry has translational symmetry in the lateral 
directions, i.e., for flat surfaces. Then the integrals in Eq. (|23|l can be easily computed, and the resulting kernel 
provides a compact account of Lifshitz theory as we will show in the next Section. Even if the surfaces are deformed 
the kernel can be obtained explicitly if one considers the limit of ideal metals, i.e., a diverging dielectric function 
e(iC). In this particular limit both pa and the operators ^'^'^ become independent of the lateral momentum ky. 
Therefore, after parameterizing the surfaces so that 2/3, j/3 are replaced by functions of the lateral coordinates yy, yy, 

respectively, the integrals in Eq. H23() correspond to Fourier transformations with respect to the lateral coordinates, 
and the kernel assumes a simple form in position space as we will demonstrate explicitly below. However, any kind 
of deviation from flat surfaces (even for ideal metals) renders A4 non-diagonal and makes the evaluation of Eq. H25|l 
a hard problem. There are basically two approaches to tackle this problem. First, one can consider the amplitude of 
the surface deformations as small compared to both the mean surface distance and other characteristic lateral length 
scales as, e.g., the roughness correlation length. Then one can apply perturbation theory to obtain the force in powers 
of the deformation profiles /iQ,(xy). This program has been carried out in detail for ideal metals in Refs. p3. l23|. 
Second, one can try to compute the force exactly by a numerical algorithm. For periodically deformed (corrugated) 
surfaces of ideal metals it has been demonstrated recently that the corresponding kernel can be transformed to a form 
which is particularly suited for an efficient numerical evaluation of the force |2j| . We expect that these techniques can 
be applied to the general case of deformed surfaces of dielectric media using our approach of expressing the force in 
terms of a kernel [Eq. (|23|) ] which contains all information about material and geometrical properties of the surfaces. 
As for ideal metals, the kernel is proportional to the vacuum Green function which, however, is now dressed by the 
operators ^'^^ which contain the reflection and refraction properties of the material. 
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III. FLAT SURFACES 



As a simple application of our approach, we consider in this section the case of flat surfaces of general dielectric 
media. In this particular limit the force between the surfaces is well-known from more conventional approaches. 
The corresponding result is known as the so-called Lifshitz theory for molecular forces 0|. In the following we will 
show that our path-integral approach provides a compact derivation of the Lifshitz result without the need to solve 
Maxwell's equations with a random source explicitly and to calculate the expectation value of the stress tensor. In 
the flat surface limit, the surfaces Sa are parameterized by (y||,i?ct) with = 0, H for a = 1, 2, respectively. Due 
to the translational symmetry of the problem, it is convenient to work in momentum space. Using the representation 



=»qiiy 



g-p(C,qii)|y3| 
MCqy) 



(26) 



of the vacuum Green function in Eq. 123() with p{(, q|| ) = i /C^ + yields 



iqiKyii-y;: 



-piC^iOlvs-y'sl 
2p(C,q||) ' 



(27) 



where we made use of the surface normal vectors n" = (0, 0,(— 1)" ^) for flat surfaces. The differential operators 
and J^"^^^ can now be expressed in momentum space with the replacements Vii ikn and Vn ih',, yielding 



-p(C,k||)|a3— yal 



MCM) 



v'a=H0 



(28) 



where we separated the factor 77"^ ^ (_i)a+/3g-[(-i)°Po(C,k||)i/o-H(-i)''p,5(C',k;|)ff^] ^j^-^^ discuss below. The 

differential operator acquires now the form 



^'"(Ck||) = 



f-k2 

ki 

V 



fclfc2 

§^[(-l)>„+93] 



: (-l)°Po 



k2\ 



fcifc2 



^,^[(-i)>a + a3] J 



ik 



(29) 



and the primed adjoint operator acts via 83 on y'^. Before we calculate from this expression the free energy and force 
between the surfaces, it is instructive to examine the structure of the above matrix. It is easily checked that the third 
row of the matrix-operator can be expressed in terms of the other two rows via {~l)"ipa^^l^ — ki^f^ + k2^2p^ ■ 
The physical reason for this can be easily understood. The relation between the rows reflects the fact that there exist 
only two independent boundary conditions for each surface. Since the surfaces are flat here, any field configurations 
can be considered as a superposition of transversal electric (TE) and magnetic (TM) modes. Each mode type 
is characterized by a scalar field which has to satisfy only one boundary condition at each surface. This will be 
demonstrated more explicitly in the Appendix where the problem is formulated from the outset in terms of two scalar 
fields representing TE and TM modes. However, for deformed surfaces this reduction to separated mode types is 
generally no longer expected to hold since the modes will mix under the scattering at deformations. 
For flat surfaces we are thus led to introduce the reduced matrix-operator 



J^«(C,k||;9,) = 



fclfc2 



fclfc2 



Ceo 

.■(-l)°pofcl 



(30) 



kl -C-^fe' + (-l)>aa,] ^ . 

which consists of two linear independent rows only. Defining J7^(C, k||) = r2"(C, ky ; ±9z) the kernel can be written 

e-p(C,k||)|2|- 



as 



.M"^'^''(C,k||;C',kji) = (2^)35(C-C')'^(')(k||-k'| ^(C, k|| ) • f2^_f'(C', kj,)) 



2p(C,k||) 



z=Hc-Hb 



(31) 
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The entries of this diagonal matrix consist of the 4x4 matrices which are given by the expression in the square 
brackets. Inserting now Eq. 13U|) into Eq. (|31|l we obtain for the expression in the square brackets the 4x4 matrix 

^iC^h) - I 1' ^""l (32) 



B A. 



2 



in terms of the symmetric 2x2 matrices 



(33) 



B 



\pI - kl-ppi){pl - kl-pp2) + kl{eie2C + kf - P1P2) ^1^2 ((ei + £2 - eie2)C^ + k| - p{pi + P2) +PiP2^ 
fcife ((ei + £2 - eie2)C^ + H ^ P(P^ + P2) 



{pI -kj- ppi){pl -kl-pp2) + kf{eie2C + k^ - P1P2) 

(34) 

using Pa = </eQ,C^ + kn , p = ./C^ + kn. The Casimir free energy per unit area can now be obtained from Eq. I|24|l . 



P1P2 



In the hmit H ^ 00 the off-diagonal elements of M in Eq. (|32|l vanish, so that we have to compute the determinant 
of the matrix 



.M(C,k||).M;,i(C,k||) = 



£1£2C' 



■BA-^ 



I 



(35) 



where / is the 2x2 identity matrix. In the above matrix we have neglected the factor rjafs appearing in Eq. (|31|l . This 
factor will have no effect on the free energy as will show at the end of this section. The determinant of the matrix of 
Eq. H35|) can be calculated using the relation 



det(y) = 1 - Tt{XiX2) + dct(XiX2) 



(36) 



for a general 4x4 matrix of the form 



Y 



' I Xi 
.X2 I , 



(37) 



Thus the free energy can be obtained by calculating the determinant of just a 2 x 2 matrix. Using Eq. H24() . the 
logarithm of the product of all the determinants for different ^„ and ky becomes a corresponding sum and integral, 
respectively, 



00 

= ^ E / ln|^(Cn,k||)^;,i(C„,k||)|, 



(38) 



where | . . . | denotes the determinant of the 4x4 matrix at fixed C, and k|| . Calculating explicitly the determinant with 
the aid of Eq. (|36|l . we obtain the final result for the Casimir or interaction free energy per unit area of the surfaces, 



-2p„g P"l -Pn Pn2 ~ Pn 
Pnl +Pn Pn2 + Pn 



-2p„H Pnl - ^iPn Pn2 ~ £2Pn 
Pnl + e.lPn Pn2 + i2Pn 



(39) 



with k = |k|||. The corresponding force per unit area is given by 



l^' f°° kdk 



Pnl + Pn Pn2 + Pn ^2p„H _ ^ 
Pnl - Pn Pn2 - Pn 



Pnl 



eiPn Pn2 + e2Pn ^2p„_ff _ ^ 



Pnl - eiPn Pn2 - ^2Pn 



(40) 



where we defined pna = \/ ^aCn + ^^'^ Pn = \/Cra + The primed sum indicates that the term for n = is to be 
multiplied by 1/2. It is important to note that the dielectric function in the above expressions is evaluated along the 
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imaginary axis only, since Cq. = ea(^C) due to the initial Wick rotation to the Euclidean field theory. Since ea(iC) is 
completely determined by the imaginary part of the dielectric function for real frequencies to, the force depends only 
the dissipative properties of the media, as expected from the fluctuations-dissipation theorem. Our result of Eqs. JSHll, 
(|40|l is in perfect agreement with the original result by Lifshitz P, Q . 

Next we consider now the zero temperature limit. This limit is obtained by the replacements Cn ^ C E'-nd 
i/PY^'n>o ~^ /o°°'^C/27r in Eqs. (|39(l and H4()|l . Following Lifshitz, we change the integration variable to 



q = Y/r+fc^7(^ and define Sa = ^g-^ — 1 + ^aiiC), yielding 
, , r°°C^dC f°^'qdq 

for the free energy, and 



-2CqH ^1 



q S2 



si+ q S2 



27r^ 



CdC / q'dq 
Ji 



si-q S2- q 



si + eiq S2 + f.2q 



gl +eig S2 + £29 ^2CqH „ ^ 

Si - eiq S2 ~ £29 



(41) 



(42) 



for the force. The above result agrees again with the original Lifshitz theory, see Eq. (2.9) in fH. 

Finally, we come back to the omitted factor t^q,^. The effect of taking into account this factor is that in the matrix- 
kernel of Eq. H32|l the matrix A2 is multiplied by e~^^^^, B is multiplied by —e^^^^ while Ai remains unchanged. 
For the matrix in Eq. H35II this means that the off-diagonal matrix ~ ^ gets multiplied by the factor — e^^^ while 
the matrix ^ A^^ gets multiplied by the inverse factor — e~^^^. Due to Eq. H36() the determinant depends only on 
the product of the two off-diagonal matrices so that the factors coming from the rjaf} drop out in the determinant of 
MM-\ 



IV. DEFORMED SURFACES OF IDEAL METALS 



In the previous section we saw that our general approach reproduces the Lifshitz theory for flat surfaces of dielectric 
media. In this section we will apply our theory to deformed surfaces. As an example we consider ideal metals with 
infinite dielectric functions • This is a reasonable approximation for surface separations which are large compared to 
the plasma wavelength of the material. However, our general result for the kernel of Eq. (|23() contains all information 
which is necessary to treat deformed surfaces of non-ideal metals or general dielectric media as well. In the latter case 
the kernel M assumes in general no particular simple form and has to be computed numerically in order to obtain the 
force. For ideal metals the kernel can be calculated explicitly and the result provides another interesting limit which 
has not been studied previously. In previous works only special deformations of ideal metals have been studied by a 
path integral approach. If the surface deformations are translational invariant in one direction as for, e.g. uni-axial 
corrugations, the electromagnetic field can be separated into TE and TM modes. This property has been used in 
[2^ I23L |23 | to describe the surface interaction by a scalar field theory. In contrast, here we will allow for general 
deformations so that no separation into TE and TM modes is possible anymore. 

Our starting point is the general result for the kernel of Eq. H23|) . After taking the limit — > 00 both pa{C,,'k^^) 

and the operators become independent of the lateral momentum kii . The kernel can then be written as 



ik||y||-ik]|y|| g-ICI [(-l)°v^i/3-(-(-l)''v^a3l 



.M"'''^"'(C,k||;C',kj|) = 27r^(C-C') / / 

X n^nf [^'="(c)-^'^^^(C')],,e(C;y-y') 



(43) 



with the differential operators 

/ ON 

= ids c 

\-id2 -C Oy 



/-ids -0 
000 

\idi CO 





f id2 





C 




—id I 


-c 







K 









(44) 



Due to the simple exponential dependence of the integrand of Eq. (|43|) on k|| it is more convenient to transform the 
kernel to position space. When we insert the height profile of the surfaces with — Ha -f ha{y\\), j/3 = -ff/3 -I- /i/5(yM ), 
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Ha = 0, -ff for a = 1, 2, the position space form of the kernel can be read off from Eq. H43() . 

X n^nf [i^'="(C)-=5f't^''(C')],,e(C;y-y') 



(45) 



Before proceeding, it is useful to discuss the two exponential factors depending on Let us start with the 

second one which depends on the height profiles ha{y\\) but is independent of the mean surface distance H. Defining 
r;Q,(C,y||) — e"!'''^"^-' v^'^atyy) the matrix has the same structure as in Eq. 1)221) but with ky replaced by y||. Due 
to the arguments given below Eq. H22|l the factors ?7Q(Ciy||) drop out of the free energy and can thus be neglected in 
the following. For the first exponential factor in Eq. H45|l this argument does not apply since the factor depends on 
H. However, we can make use of the fact that the factor does not depend on the lateral coordinates yy . The effect 
of this exponential factor is that every 2x2 sub-matrix of M resulting from keeping (j, I), (C,y||) and (C',y||) fixed 

is multiplied by the same only C, dependent factors. The two diagonal elements are multiplied by 1 and e^^'^l^^-^, 
respectively, while the off-diagonal elements are multiplied by e"'"''^^^. It is easy to check that this leads to the 
simple factor e~^l'»lv^^' for the determinant of if denotes the dimension of the matrix M. However, this factor 
will drop out when taking the determinant of the ratio . Therefore, the first exponential factor in Eq. H45|) 

can be omitted as well. 

Now the kernel assumes a simple form. Expressing the surface normal vectors in terms of the height profile by 




(46) 



with /iq, J- — djha and 5a = 1 + (Vy/ic)^ the kernel can be written as a functional of the height profile. For the same 
reason as the factors ?7Q(C,y||) could be omitted above, we can neglect the normahzation factor (— l)"/v^ of the 
normal vector. Thus we obtain for the differential operator 



[fi"^'"](C,yii) = 



iha,id3 + idi 




(47) 



We observe that the third row of the matrix in H47|l is linearly dependent since ha^i^i. 



Jlg^j. Therefore, 



in the Lifshitz theory limit discussed earlier, the matrix has to be reduced to its first two rows. The linear dependence 
of rows reflects the fact that for ideal metals there are, in fact, only two independent boundary conditions for each 
surface. As mentioned earlier, for general deformations a reduction to TE and TM modes as it appears in Lifshitz 
theory is not possible. However, for ideal metals the boundary conditions can be simply written as 



[na(y) X E(C,y) 



= 0. 



(48) 



This boundary condition requires the two tangential components of the electric field at the surface to vanish locally. 
Saying it differently, the limit of infinite conductivity converts the three originally non-local boundary condition into 
two local conditions. The final result for the matrix kernel is now given by a 4 x 4-matrix with j, Z = 1, 2, 



^"'''^■'(C,y|i;C',yji) = 27rJ(C-C')[f^"-f^"'1,,^(C;y~y') 



with the deformation dependent differential operator 



3-1 



(49) 



(50) 



acting on the vacuum Green function. Here the prime on h indicates the dependence on the primed variable yy . 
This kernel together with the formula of Eq. (|24|l yields the exact free energy of the interacting surfaces. Generally it 
is not possible to give a closed analytical expression for the determinant of M of H49|) . However, either perturbative 
[2^ I23 or numerical |2^ techniques can be used to evaluate the free energy and force from the kernel using Eq. H24|) 
and Eq. 
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V. DISCUSSION AND OUTLOOK 

We derived from a path-integral quantization of the electromagnetic gauge field with non-local boundary conditions 
an effective action for the molecular forces between dielectric media with deformed surfaces. From this effective action 
we derived the Lifshitz formula for flat surfaces and an explicit expression for the matrix kernel which determines 
the interaction between ideal metals with arbitrary deformations. We believe that our approach will be a useful 
starting point for the calculation of correlations between modifications of the ideal Casimir force of Eq. due to 
material properties (finite conductivity) and geometry. Our theory is certainly not a final answer to the problem of 
correlation effects. Rather it is intended for a generalization of the Lifshitz formula to deformed surfaces. As for the 
Lifshitz formula, for an explicit calculation of the force one, of course, has to specify a suitable dielectric function 
and, in addition, a height profile describing the geometry. An precise result for the force can presumably only be 
obtained from a numerical evaluation of the effective action. However, certain limiting cases should be accessible to 
a perturbative analysis. For example, it would be interesting to consider the cases of small and large distances H 
between the surfaces separately. For small H, the interaction is dominated by high frequencies and is sensitive to the 
model from which the frequency dependence of the dielectric function e(a;) is obtained. The correct choice of e(a;) 
for real metals is currently being discussed but this is not the matter of the present discussion. For large H a useful 
approximation is to assume e to be constant. In a large e expansion a strong geometry dependence and non-additivity 
of Casimir forces between ideal metals should be found. In the opposite limit of rarefied media the non-additivity 
effects and thus the sensitivity to geometry should decrease if e approaches one. Another point of recent dispute is 
the choice of boundary conditions for real metals |23|. Although this is not the subject of the present paper, we note 
that our path integral approach is sufficiently general to include any type of boundary conditions as, e.g., the so-called 
impedance boundary condition ||4^ ^] . 

APPENDIX A: SCALAR FIELD THEORY FOR FLAT SURFACES 

The scope of application of scalar field theories is more restricted than that of the gauge field formalism, but 
calculations become simpler in many cases. An interesting example is the case of flat surfaces of dielectric media. In 
Section ITTll we re-derived the Lifshitz theory for dielectric media. There we argued that the field can be decomposed 
into TE and TM modes each of which have to obey only one boundary condition on each surface, leading to a reduced 
matrix kernel. In this Appendix we show explicitly that the Lifshitz theory can be obtained directly from a path 
integral quantization of a scalar field which fulfills a suitable boundary condition. Therefore, no gauge fixing procedure 
is needed in this special situation. 

In order to decompose the electromagnetic field into transversal electric (TE) and magnetic (TM) modes, we have to 
specify a preferred spatial direction. Due to the rotational symmetry in the lateral plane of the flat surface geometry, 
we can choose without any restrictions the j/2 direction. Actually, it turns out that it is useful to adapt the choice of 
the direction to the lateral momentum ky of the field mode. Note that we can do this since the modes for different 
k|| are decoupled for fiat surfaces. Therefore, in the following we will choose for a given ky the lateral coordinates 
so that ^2 = 0, and y2 defines the longitudinal direction. Following the parameterization of the TE and TM modes 
for waveguides |4^, the longitudinal components of the electric (TM modes) and magnetic (TE modes) field define a 
(real- valued) scalar field $, 

E2 = ^ B2=Q for TM modes 

(Al) 

^2 = $ E2 = for TE modes. 
The transversal components Ej = Ej{(^, k|| 

El 

Ez 

for TM modes, and 

El 
E3 

for TE modes, where we performed already a Wick rotation to imaginary frequency, uj — > Using Maxwell equations 
it can be shown that the dynamics of the scalar field <f> are governed by the usual wave equation, corresponding to 



z) and Bj = Bj{(, k|| ; z) of the electromagnetic field are then given by 
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the Euclidean action 



5e[1>] = (V$)2 (A4) 

with the partition function 

Zo = / f^$e-^^[*l. (A5) 



In this section we consider only the zero temperature case. Finite temperatures can be treated in analogy to the gauge 
field approach in section FlI Bl by introducing Matsubara frequencies. The boundary condition for $ can be derived 
from the general condition wc found in section FlI Al From Eq. pi|l we obtain, using the Fourier representation of the 
material Green function Q'^ , in Euclidean space the conditions 

-C^o.{^C,)E2-^klB■i-{-lrp^B^ = (A6) 

Cea(«C)^^i-*feS3-(-l)>aB2 = (A7) 

ik|| • B|| - (-1)>„S3 = (A8) 



for the flat surface Sa for a = 1,2, and with pa = y^ea(iC)C^ + k|. Now we make use of the fact that we can constrict 

the analysis to the case ^2 = by a suitable choice of the lateral coordinates. After inserting the electromagnetic field 
components as given by Eqs. (IA1|) . IIA2p and ljA3|l . the above conditions collapse to a single boundary condition for 
the scalar field $. Depending on the type of mode, we obtain the condition 



1 - r„ a,. 



$U=ff„ = (A9) 



with 



r„ = — for TM modes, r„ = — for TE modes, (AlO) 

Pa Pa 

where = {—l)°'^^dz denotes the normal derivatives of the surfaces. Due to the ky -dependence, the condition 
(|A9|I is non-local in position space. In the limit of ideal metals, — > co, the above condition reduces to the well 
known Dirichlet and Neumann boundary conditions for TM and TE modes, respectively. We note that Fq, is real, 
since the same is valid for Eq, on the imaginary frequency axis. Therefore we need to consider only a real-valued field 
$. The restricted partition function for this field reads for both types of modes 

Z = 2o-' /'i^<i>e-^-Wnn'5[(l-r„^„J$U=^J, (All) 

" C,k|| 

where we implemented the boundary constraints again by delta- functions. Now we proceed in analogy to the treatment 
of the gauge field path integral in section III Bl We introduce two auxiliary field -ipa, one for each surface, in order to 
replace the delta-function. After integrating out <&, we obtain for the partition function 



with the effective action 



Sosm^l [ I V'a(C,k||)X"'^(C,k||;C',k|j)#(C',kj|) (A13) 



and the 2x2 matrix kernel 



^ " ^ " 2p l(l_pri)(l-pF2)e-P^ l~p^Vl I 



with k ~ |k|||, p — y^Q^ + k^ and a summation over a and f3 is implicit. From this kernel, the Casimir energy can 
be quite easily calculated. Using Eq. H24f) in the zero temperature limit, we obtain for the Casimir energy per unit 
surface area 



l+pFi l+pFa 



(A15) 
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where | . . . | denotes the determinant of the 2x2 matrix at fixed C and ky. The matrix Moo is diagonal since for 
if — > oo the off-diagonal elements in Eq. (|A14I) vanish. Therefore, the determinant of the matrix product can be 
easily obtained, yielding the last expression in Eq. I)A15|I . If we substitute Fq, according to Eq. IjAlOjl in Eq. I|A15|I 
by 1/pa or by tajVa we obtain the contribution of the TM or TE modes, respectively, to the energy. It is easily seen 
that the sum of the energies from both modes reproduces the Lifshitz theory result we obtained before in section UlTI 
from the gauge field approach, cf. Eq. (jSHJ- 
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